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Abstract 

We find an analytic solution of the Bethe Ansatz equations (BAE) for the special 
case of a finite XXZ spin chain with free boundary conditions and with a complex 
surface field which provides for U q (sl(2)) symmetry of the Hamiltonian. More precisely, 
we find one nontrivial solution, corresponding to the ground state of the system with 
anisotropy parameter A = \ corresponding to = — 1. 

With a view to establishing an exact representation of the ground state of the finite 
size XXZ spin chain in terms of elementary functions, we concentrate on the crossing- 
parameter n dependence around n = tt/3 for which there is a known solution. The 
approach taken involves the use of a physical solution Q of Baxter's t-Q equation, 
corresponding to the ground state, as well as a non-physical solution P of the same 
equation. The calculation of P and then of the ground state derivative is covered. 
Possible applications of this derivative to the theory of percolation have yet to be 
investigated. 

As far as the finite XXZ spin chain with periodic boundary conditions is concerned, 
we find a similar solution for an assymetric case which corresponds to the 6-vertex 
model with a special magnetic field. For this case we find the analytic value of the 
"magnetic moment" of the system in the corresponding state. 

Dedicated to Rodney Baxter 
on the occasion of his 60th birthday. 
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I. Introduction 



It is widely accepted that the Bethe Ansatz equations for an integrable quantum spin chain 
can be solved analytically only in the thermodynamic limit or for a small number of spin 
waves or short chains. In this paper, however, we have managed to find a special solution of 
the BAE for a spin chain of arbitrary length N with N/2 spin waves. 

It is well known (see, for example |L] and references therein) that there is a correspondence 
between the Q-state Potts Models and the Ice- Type Models with anisotropy parameter 
A = y y-. The coincidence in the spectrum of an N-site self-dual Q-state quantum Potts 
chain with free ends with a part of the spectrum of the U q (sl(2)) symmetrical 2N-site XXZ 
Hamiltonian ([]]) is to some extent a manifestation of this correspondence. 

H xxz = Wn^n+l + °n °"n+l + °"n°"n+l + ~ T (°"n ~ °"n+l)}> CO 

where A = (q + q~ 1 )/2. This Hamiltonian was considered by Alcaraz et al. [0 and its 
U q (sl(2)) symmetry was described by Pasqier and Saleur ||. The family of commuting 
transfer-matrices that commute with H xxz was constructed by Sklyanin || incorporating a 
method of Cherednik 0. 

Baxter's T-Q equation for the case under consideration can be written as || 

t{u)Q{u) = <f>(u + r]/2)Q{u -7]) + cj)(u - r]/2)Q{u + rj) (2) 

where q = expir), (f>(u) = sin 2u sm 2N u and t(u) = sm2uT(u). The Q(u) are eigenvalues of 
Baxter's auxilary matrix Q(u), where Q(u) commutes with the transfer matrix T(u). The 
eigenvalue Q(u) corresponding to an eigenvector with M = N/2 — S z reversed spins has the 
form 

M 

Q{u) = J\ sin(w - u m ) sin(w + u m ). (3) 

m=l 

Equation (||) is equivalent to the Bethe Ansatz equations 

sin(n fc + fj/2) 
_sin(M fc - rj/2) 

Baxter's equation can be interpreted as a discrete version of a second order differential 
equation|7|, ||. So we can look for its second independent solution P{u) with the same 
eigenvalue T(u). 

In a recent article U Belavin and Stroganov argued that the criteria for the above men- 
tioned correspondence is the existence of a second trigonometric solution for Baxter's T-Q 
equation and it was shown that in the case rj = 7r/4 the spectrum of H xxz contains the 
spectrum of the Ising model. In this article we limit ourselves to the case 77 = ir/3. This case 
is in some sense trivial since for this value of r], H xxz corresponds to the 1-state Potts Model. 
We find only one eigenvalue Tq(u) of the transfer-matrices T(u) when Baxter's equation (0) 
has two independent trigonometric solutions. Solving for T(u) = Tq(u) analytically we find 
a trigonometric polynomial Qo(u) the zeros of which satisfy the Bethe Ansatz equations 



M 

n 



u, 



+ rj) sm(u k + U m + 7]) 



sin(u fc -u m -rj) sm(u k + u m - 77) ' 



(4) 
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(0). The number of spin waves is equal to M — N/2. The corresponding eigenstate is 
the groundstate of H xxz with eigenvalue E Q — |(1 — N), as discovered by Alcaraz et al. 
numerically. 

When the first version of this paper 



15] was reported at the Workshop in ANU (February 



2000), Baxter showed us his article [TIJ where he writes about the existence of a simple 
eigenvalue of the T-matrix for the 8- vertex model for the special case /i = tt/3. About 30 
years ago Baxter discovered [17| that the ground state energy of the XYZ-spin chain which 
is described by the Hamiltonian 



N 



H 



xyz 



(5) 



ra=l 



has an especially simple value E = 
J z satisfy the constraint J x J y + J y Jz 



when J x , J y and 



(J x + J y + J z )/4: for the case jj, = 7r/3 
JzJx = 0. This value was found in the thermodynamic 
limit only. Later in 1989 |16[ Baxter considered a determinant functional relation and found 
a very simple solution for eigenvalue of T-matrix (see Section II below) for the case \i = 7r/3. 

There are however problems with the realization of this simple solution. If we consider, for 
example, the usual XXZ-spin chain with periodic boundary conditions we will find, solving 
Baxter's T-Q equation, that the degree of Q(u) is N/2 + 1 (N even) and there is no Bethe 
vector for this simple solution. We believe that for the case without a field, a "physical" 
solution Q(u) must have degree less than or equal to N/2. Baxter informed us that the 
addition of a special field saves the simple solution and below we find the corresponding 
Q(u). 

At first glance it seems we have a clear picture. For the open chain with special boundary 
conditions we have central charge c = (see for example |TB[), so finite-size corrections are 
absent and we have a simple solution for the ground state. On the other hand, for periodic 



boundary conditions without a field we have due to Hamer |19j c ^ 0, so there is no simple 
solution. 

We wish to stress that essentially only chains with even length were considered. Hamer |19 
considered only even N. However, for odd N, one can find []20| a simple eigenvalue To = 
(a + b) N for the case when the weights of the 8- vertex model satisfy the condition 



(a 2 + ab)(b 2 + ah) = (c 2 + ab){d 2 + ab), 



(6) 



which is equivalent the condition \i = %/3. 



II. Consequences of the existence of the second "inde- 
pendent" periodic solution 

This question was considered in article |9| . Here we use a variation more convenient for our 



goal. This section is mainly due to Baxter [|1(| (see also ||21|| ). 

Let us consider T-Q equation (0) for rj = ^, where L > 3 is an integer. Let us fix a 
sequence of spectral parameter values Vk = v + r/k, where k are integers and write (pk = 
4>{ v k — ^/2), Qk = Q( v k) and tk = t(vk). Functions 0(w), Q{u) and t{u) are periodic 
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with period tt. So the sequences we have introduced are also periodic with period L, i.e., 
4> k+L = fc , etc.. 

Setting u = v k in (|2]) gives the linear system 



tkQk — 4>k+lQk-l + 4>kQk+l- 

The matrix of coefficients for this system has the three-diagonal form 



(7) 





-to 


00 


.. 
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01 




02 


~tl 


01 .. 


. 





M = 





03 


-t 2 ■■ 


. 







0L-1 





.. 


• 00 


-t L - 



Taking ^ where m is an integer, we have 0& 7^ for all k. 

It is straightforward to calculate the determinant of the L — 2 x L — 2 minor obtained 
by deleting the two left most columns and two lower most rows. It is equal to the product 
—0^0203 . . . 0l-i, which is nonzero, hence the rank of M cannot be less than L — 2. Here 
we are interested in the case when the rank of M is precisely L — 2 and we have two linearly 
independent solutions for equation (|7|). Let us consider the three simplest cases L = 3, 4 
and 5. Parameter r\ is equal to |, ~ and | respectively. 

For L = 3 the rank of M is unity and we immediately get to = — 02, t\ = — <po and 
t>i = —0i- Returning to the functional form, we can write 



T (u) = to(u)/ sin2w = — <p(u + — )/ sin2w 



2N 
COS u. 



(9) 



This is the unique eigenvalue of the transfer-matrix for which the T-Q equation has two 
"independent" periodic solutions. It is well known (see, for example, 0) that the eigenvalues 
of H xxz are related to the eigenvalues T(u) by 



E 



T'0n/2) 

N cosr] + (— — -— + tan 77) sin 77. 



(10) 



Tfo/2) 

For the eigenstate corresponding to eigenvalue @ we obtain 

E = l(l-N). (11) 

This is the groundstate energy which was discovered by Alcaraz et al. numerically. 

In the next section we find all solutions of Baxter's T-Q equation corresponding to T(u) = 
Tq{u). In particular we find Q(u), the zeros of which satisfy the BAE ([|). 

For T] = I and transfer- matrix eigenvalue Tq(u) = cos 2N u, T-Q equation (g) reduces to 



<j)(u + 3r}/2)Q{u) + (j>{u - v/2)Q{u + 77) + 0(w + r]/2)Q{u — 77) = 0. 



Keeping t(u) arbitrary, for the moment, we can rewrite 
t(u) = 



as 



u + W 2 ) — 7TT~\ ^ 0( M _ W 2 )" 



Q{u) 



Q{u) 



(12) 



(13) 
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We also have 



t(u + r)) = 4>(u + 3r]/2)- 



Q(u) 



+ 4>(u + 77/2) 



Q(u + 2r]) 



Q{u + rj) rv Q(u + v) ' 

Multiplying these equations together we obtain the fusion relation 

t{u)t{u + r]) = 4>(u - r)/2)(f)(u + 377/2) + 4>(u + r]/2)i{u), 

where 

Q(u-rj) 



(14) 



(15) 



t(u) 



Q{u)Q{u + 77 



{<j>{u + 3r]/2)Q(u) + <P(u - r]/2)Q(u + rj) + </>(u + rj/2)Q(u - 77)). (16) 



In the case under consideration we have t(u) = 0. Fusion relation ( fTS] ) reduces to the simple 
equality 

t (u)t (u + r]) = 4>{u - 77/2)0(m 



3^7/2), 



(17) 



which is a kind of inversion relation. 
For L = 4 



M 



-to 


00 





01 
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~tl 


0i 








03 


-t 2 


02 


03 





00 





Deleting the second row and the forth column we obtain a minor with determinant — 0o03(^o+ 
£2). It is zero when t 2 = —to, i.e., t(u + f ) = —t(u). Considering the other minors we obtain 
the functional equation 



t(u + n/8)t(u - tt/8) = 0(tj + 7t/4)0(tj - tt/4) - 0(7j)0(7J + tt/2). 



(19) 



This functional equation was used in to find t(u) and show that this part of the spectrum 
of H xxz coincides with the Ising model. It would be interesting to find a corresponding Q(u). 

Lastly for L = 5, minor M 35 (the third row and the fifth column are deleted) has deter- 
minant 0004(^1 + 01^3 — 0002)- Setting this to zero we have 



t(u)t(u + tt/5) + 0(u + 7r/10)t(u + 3tt/5) - 0(n - tt/1O)0(t7 + 3tt/10) = 0. 



(20) 



It is not difficult to check that in this case all 4 x 4 minors have zero determinant and that 
the rank of M is 3. Thus we have two "independent" periodic solutions of Baxter's T-Q 
equation. 

Note that this functional relation ( |2"0"| ) coincides with the Baxter-Pearce relation for the 
hard hexagon model. The connection between ( f20|) and a special value of the rank of the 
matrix of coefficients for system ([7p was remarked upon in ]n]] by Andrews, Baxter and 
Forrester. 

For general L we obtain the same truncated functional relations that have been obtained 
in with the same assumptions. Note that for the ABF models which are a general- 
ization of the hard hexagon model, the truncated functional relations have been proved by 
Behrend, Pearce and O'Brien in fT^fl . 
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III. Solution of Baxter's Equation for 77 = | and T = Tq 

Equation ( |T2"D can be rewritten as follows: 

lir Ait 

/(„) + /(« + T ) + /(« + T ) = 0, (21) 

where /(f) = sinf cos 2Ar (f/2) Q(v/2) has period 27r. The trigonometric polynomial /(f) is 
an odd function so it can be written 

K 

/(f) = E c fc sin ^ ( 22 ) 
fc=i 

where K is the degree of /(f). Equation (|21~D is equivalent to 

c 3m = 0, me Z. (23) 
The point f = 7r is a zero of /(f) of order 2N + 1, so we obtain 

(-f)7(«)Ur = 0, z = 0,l,...,2iV. (24) 

Of 

For even 2 this condition is immediate, whereas for i = 2j — 1 we use (^) to obtain 

£ (-l) k c k k 2 ^ = 0, j = 1,2,..., iV. (25) 

fc=l,fc^3m 

Our problem is a special case of a more general problem which can be formulated as 
follows. Given a set of different complex numbers X = {xi,x 2 , ■ ■ ■ ,xi} we seek another 
complex set B = {f3\, /3 2 , ■ ■ ■ , fti} where ft ^ for some i, so that 

i 

52PiP(xi)=0 (26) 
i=i 

for any polynomial P(x) of degree not more than N — 1. It is clear that for I < N the 
system B does not exist. If j3i ^ 0, for example, the product (x — x 2 )(x — x%) . . . (x — xj) 
provides a counterexample. 

Let I — N + 1. We try the polynomials 

N 

Pr= II (x-Xi), r = 1,2,..., AT. (27) 
Condition (|26"D gives [3 r P r {x r ) + PjP r (xj) = and we immediately obtain 

N+l 

Pr = COnst Yl ( X r — ^i) 1 , (28) 
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which is a solution because the system ( p7|) forms a basis of the linear space of TV — 1 degree 
polynomials. So for / = N + 1 we have a unique solution (up to an arbitrary nonzero 
constant) given by (|28|). 

Returning to ( ^5|) we first consider N = 2n, n a positive integer. Fix I — N+l — 2n + l. 
The degree K becomes 3n + 1. It is convenient to use a new index k = \3k + 1|, where 
\k\ < n. Equation ( p5|) can be rewritten as 



£ &(3« + l) 2 Ci-U = 0j J = 1, 2, . . . , A^, (29) 
where we use new unknowns (3 K = (— 1) k C| 3k+ i| |3/c + 1| instead of c^. Using (p8|) we obtain 

n 

= const J] ((3/« + l) 2 -(3p+l) 2 )- 1 . (30) 

We can rewrite this using binomial coefficients as 

ft = const' ( K+ i)( 2 ;_ + j)( 2 ; + -j). (31) 
The old system of unknowns is given (up to an arbitrary constant) by 

c M = (-ir S gn ( K + i)( 2 ;_ + j)( 2 ;;j) (32 ) 

and using (F22|) we obtain the function /(v), 



/(.) = tj-lT ( 2 „"_ + j ) ( l" + j ) + D«- (33) 

We recall that the solution of Baxter's T-Q equation for T{u) = T (u) is given by 



QoM= ■ / (2 "L (3D 
sin 2m cos /iV it 



and its zeros {m^} satisfy the BAE 



IV. Distribution of zeros of Qo(u] 

Function f(v) p3] ) can be shown to satisfy a simple second order linear differential equation. 
The coefficient functions of this equation are closely connected with the density function of 
the zeros of Qq(u) in the thermodynamic limit. 

Let us introduce x = exp (3iv) and rewrite f{y) as 

n 

f(v) = F(x) = F+(x) - F~(x), F + (x) = V a K x K+ i F~(x) = F + (l/x). (35) 
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We have a multiplicative recursion relation for a K 

a K+ i (n — k){ti — k — 2/3) 



(36) 



a K (n + K + l)(n + k + 5/3) 

which gives [] 

n 

Y W+i( n + K + l)(n + k + 5/3) + a K (n - n)(n - k - 2/3)}x K+ ^ = (37) 



or 



it it 

Y a K {(n + k + 1/3) 2 - l/9}x K ~l + ^ a K {(n - k - 1/3) 2 - l/9}x K+ ^ = 0. (38) 

K=—n K=—n 

Using the operator 9 = x-£; we can rewrite this as 

{((6 + nf - l/9)/x + (0 - n) 2 - 1/9}F+ = 0. (39) 

This equation is invariant under the transformation x — > 1/x, so F~(x) = F + (l/x) is also a 
solution for (j3~9| ). 

Alternatively, we know that in a neighbourhood around the singular point x — — 1 there 
are two solutions of ( p9|) which behave as (x + l) a (l + 0(x + 1) + . . .). Substituting this into 
(^9|), we obtain a.i — and «2 = 4n + l. This method can be used as an another way to find 

/(*)• 

Returning to variable v we have 

{exp{-3iv)({d/dv + 3m) 2 + 1) + ((d/dv - 3m) 2 + l)}f{v) = 0. (40) 
Multiplication with exp(3w>/2)/ cos(3f/2) gives 

d 2 f df 

-4 + 6ntan(3v/2)-^- + (1 - 9n 2 ) f = 0. (41) 
dv 1 dv 

The zeros of f(v), the density of which is important in the thermodynamic limit, are 
located on the imaginary axis in the complex f-plane. So it is convenient to make the change 
of variable v = is. It is also useful to introduce another function g(s) = f(is)/ cosh 2n (3s/2). 
The differential equation for g(s) is then 

y V 2 cosh 2 (3s/2) / y v ; 

Let g(so) = 0. For large n we have in a small vicinity of sq 



«" I J" ,J 2 g = 0. (43) 



cosh(3s /2) J 



1 F + (k) = F(— 2n, 2/3— 2n, 5/3, — x)x 1 ^ 3 n , where F(a, b, c, x) is the usual Gauss hypergeometric function. 
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This equation describes a harmonic oscillator with frequency uo = 3n/ cosh(3so/2). The 
distance between nearest zeros is approximately As = tc/uj and we obtain the following 
density function which describes the number of zeros per unit length 

p(s) = 1/As = uj/tt = — — — . (44) 

HK ' ' 1 vrcosh(3s/2) v ; 

It is possible to use the theory of Sturm-Liouville operators for a more careful consider- 
ation of (|42|). We note that this equation has a history as rich as the BAE. Eckart [13|] used 



the Schrodinger equation with a bell-shaped potential V(r) = —Gj cosh 2 r for phenomeno- 
logical studies in atomic and molecular physics. Later it was used in chemistry, biophysics 
and astrophysics, just to name a few. For more recent references see for example [ p}] . 

V. ^-dependence of the ground state energy of the finite 
XXZ spin chain 

We recall Baxter's t-Q equation in terms of the second solution P, 

t(u)P{u) = <f)(u + r)/2)P(u -rj) + <t>{u - r,/2)P(u + rj) (45) 

P(u) is characterised by having degree N + 2. As in the case of Q{u) we set up a function 
f{y) = smv cos 2iY (t>/2) P(v/2), and write 

K 

f(v) = ^2c k smkv, (46) 

k=l 

where K is now 3n + 2. The number of equations for coefficients Ck is determined by t{u) 
and is iV as before. Thus in this case we have freedom in two of the c^. 

Solving equation (p5|) with this extra freedom, we obtain f{y) = afp(v) + /3/g(f ) where 
a, (3 are arbitrary, fqiv) is the solution for Q(u) given by ( |33D and 

/,(«) = ± (-!)■( 2 n "_ + j ) ( 2 ; + " j ) S in( 3K + 2), (47) 

K=—n x / \ < / 

Baxter's equation can be interpreted as a discrete version of a second order differential 
equation, so we can express its coefficients in terms of the two independent solutions JF|, 

(j){u) = P(u + r]/2)Q(u - Tj/2) - P(u - r]/2)Q(u + rj/2) = p+q_ - p_q + (48) 

and 

t{u) = P(u + r])Q(u -T))- P(u - 7])Q{u + rj) = P+Q_ - P_Q+. (49) 
Consider now the functions t, Q and P as functions of two variables u and rj. <p{u) does 



not depend on rj, so taking the derivative w.r.t. rj from (HBI) we obtain 



= 5p + q_ + p + 5q_ - 5p_q+ - P-5q + + ^{p'+Q- ~ P+Q- + P-Q+ ~ P-S+)- ( 50 ) 
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We use 5 for the derivative w.r.t. 77 and / for the derivative w.r.t. u. Equation ( pE9|) leads to 

St = 5P + Q_ + P + 5Q_ - 5P_Q + - P_<5Q + + P|Q_ - P + Q'_ + P'_Q+ - P_Q' + (51) 

For 77 = Ti / 3 we have, for example, p + = P(w + 7r/6) = P{u + n/2 — n/3), so if we shift the 
variable u in (^) by 7r/2 we get 



= 5p_q+ + p_5g + - 5p + q_ - p + 5g_ + -(pLq+ - p+g'_ + p;<2- - p-Q+) (52) 



1 

2 1 

Adding this equation to (^1]) we arrive at the formula for the derivative of the largest eigen- 
value of the transfer matrix 



7]=w/3 



-(P;Q_ - P + Q'_ + P'_Q + - P_Q' + ) (53) 



One can use ( p3|) and fl47j) to get an explicit expression. 

Note that "Wronskian" (^) reduces to a complicated combinatorial identity 

*-* - 1) - - 1) - (-D^r 2 " : 1 ) ( 2 " : 1 ) ( J n l 2 , ), (54) 



2n + 1 V 2n / V 2n J \ K + 2n + 1 
where K is an integer 1 — 2n < K <2n — 1 and 



" n ' 2n + | \ / 2n - f \ / 2n + f \ / 2n - f 



k=— n K=—n 



VI. Finite XXZ spin chain with periodic boundary con- 
ditions and with additional magnetic field 

The assymetric 6- vertex model under consideration was investigated by Baxter [p2 ]. The 
associated spin chain Hamiltonian can be found in the paper of Perk and Schultz |[23|| . 
Baxter's T-Q equation with additional field can be written as follows 

T(u) Q{u) = exp(h) sm N (u — 77/2) Q(u + 77) + exp(-/V) sin 7V (u + r//2) Q(w - 77). (56) 

Using a similar procedure to Section II (see also Baxter's original paper |TB[) we find that for 
77 = 27r/3 or for 77 = tt/3 we again have a simple solution for a special choice of the field. Let 
us fix 77 = 2-7r/3 and h = it]. For this case the transfer matrix T has an eigenvalue — sin n u. 
Using the method of Section III, we find that 

f n — i \ 1 n 

K=0 



K=0 

where n = N/2. 



f(u) ee sin N u Q(u) = \ 1 ) [ n _ \ ) «p(u*(6k - 3n)) 

-j:( n ~J)( n n _~ 1 l K ) eM*u(-6n + 3n- 2)) (57) 
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The trigonometric polynomial P(u) = Q(—u) satisfies a conjugate equation which can 
be obtained from ([56]) by changing the sign of parameter h. We easily obtain (see for 
example M) the analogues of (EST) and (E3): 



(j){u) = exp(-h)P(u + r)/2)Q{u - rj/2) - exp(h)P(u - r)/2)Q(u + r)/2), (58) 

and 

T{u) = exp(-2h)P(u + rf)Q{u -rj)- exp(2/i)P(u - r])Q(u + rf). (59) 

Using similar procedure to the previous section we find the derivative of eigenvalue T(u) 
w.r.t. parameter h, which is the "magnetic moment" we refered to in the abstract. 



dT 
~dh 



= -3{exp(/i) Q{-u - 7])Q{u -v)+ exp(-/i) Q(-u + rf)Q{u + rf)}. (60) 

fe=2i7r/3 
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